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0\ ; Abstract 

C ■ If the zero-field transition in high temperature superconductors such as 

^ . YBa2Cu307_5 is a critical point in the universality class of the 3-dimensional 

y_^ , XY model, then the general theory of critical phenomena predicts the ex- 

^^1 istence of a critical region in which thermodynamic functions have a char- 

c/5 ' acteristic scaling form. We report the first attempt to calculate the univer- 



C^ 



sal scaling function associated with the specific heat, for which experimental 
data have become available in recent years. Scaling behaviour is extracted 
from a renormalization-group analysis, and the 1/A^ expansion is adopted as 
^ \ a means of approximation. The estimated scaling function is qualitatively 

O \ similar to that observed experimentally, and also to the lowest-Landau-level 

scaling function used by some authors to provide an alternative interpretation 
of the same data. Unfortunately, the 1/iV expansion is not sufficiently reliable 
^ , at small values of A^ for a quantitative fit to be feasible. 
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I. INTRODUCTION 

In recent years, a considerable body of experimental evidence has accumulated to suggest 
that the zero-field transition in certain high-temperature superconductors, most notably 
YBa2Cu307_5 (YBCO), is a critical point in the universality class of the 3-dimensional XY 
model [0-0] • If this is the case, then, in the presence of a sufficiently small magnetic field B, 
the specific heat is expected to have a singular part which exhibits the scaling behaviour 

aing(T,i?) = i?-"/2T(x), (1) 

where a ~ —0.013 and u ~ 0.67 are critical exponents and the scaling variable is x = (T — 
Tc)B~^^'^'^. Similar scaling forms are expected for other thermodynamic quantities, such as 
the magnetization. In the limit S — > 0, the scaling function must behave as C{x) -^ C±\x\~"', 
so that Csing(T, 0) = C±\T- T^l-"", where + refers to T > T^ and - to T < T^. For YBCO, 
zero-field measurements of the specific heat presented by several authors seem to agree 
well with this prediction m| and to be consistent with the universal values of a and of the 
amplitude ratio C+/C- as determined by precision measurements of the superfiuid transition 
in "^He, which is also in the universality class of the 3-dimensional XY model |]9|Jl0[]. A claim 
has recently been made that the zero field singularity is actually characterized by different 
exponents a;+ and a;_, above and below T^., which would not be consistent with any ordinary 
type of critical point [|rT|]. It has also been argued, though, that this conclusion rests on an 



inappropriate background subtraction |12 



In a nonzero applied field, one can test the scaling form (|I|) by the extent to which data 
for i?"/^'^Csing(T, B) collapse to a common curve when plotted as a function of x. Here, 
matters are complicated by the fact that a different kind of scaling behaviour 

C(r,5)^CLLL(xLLL) (2) 

is expected when only the lowest Landau level is significantly occupied [|l^,|l^]. Here, the 
scaling variable is xlll = {T — Tc2{B)) / {T B)"^^^ , where T — Tc2{B) or, equivalently, B = 
Bc2{T) is the upper critical field of the Ginzburg-Landau theory. Since a in (|I]) is very 
small, and l/2z/ ^ 0.75, the two predictions are rather hard to distinguish. Some authors 
claim that lowest-Landau-level scaling works just as well as, or indeed better than, critical 
point scaling [|l^-|2l[. For HgBa2Ca2Cu308+5 [|^, specific heat data appear to collapse to a 



common curve when plotted in the form of ([^), but the scaling function C{x), which ought 
to be universal, is apparently rather different from that found for YBCO. For LuBa2Cu307, 
the authors of Ref. [^ find that a two-dimensional lowest-Landau-level scaling form best 



fits the data, though they claim that it is also consistent with 3-dimensional XY scaling for 
fields below IT. Most recently, Junod et al [^ have concluded that optimally-doped YBCO 
is the only material to show convincing evidence of critical-point scaling. 

Theoretically, it seems that the scaling form (^ is an unambiguous prediction of the 



theroy of critical phenomena [^ and ought to be observed sufficiently close to the zero-field 
critical point. Lowest-Landau-level scaling, on the other hand, is to be expected in large 
fields, in the neighbourhood of the upper critical field. There is in principle no region where 
both scaling forms could be simultaneously valid [^. There is, however, no reliable means 



of estimating the largest field in which critical-point scaling ought to be observable or the 



smallest field consistent with lowest-Landau-level scaling. Calculations are somewhat sim- 
plified by the lowest-Landau-level approximation, and scaling functions have been obtained 
by both perturbative p7|,p9|,p8| and nonperturbative BO] methods. In particular, Tesanovic 



and Andreev ^^ have obtained scaling functions which agree quite well with experimental 



data for the specific heat and magnetization of YBCO, though the fit is rather better in the 



case of the magnetization than the specific heat [|T8|. 

For the critical-point scaling function, no theoretical estimate has been obtained (al- 
though some general consequences of scaling have been discussed in [0), and the calcu- 
lation of this scaling function is the object of the work reported here. The calculation is 
based on the Ginzburg-Landau- Wilson model of an isotropic superconductor. Although the 
superconductors of interest are anisotropic, layered materials, this seems to be a reasonable 
approximation in the case of YBCO. More generally, in fact, it is the divergence of the 
coherence length near a critical point which gives rise to characteristic critical phenomena. 
To the extent that the critical behaviour is that of a 3-dimensional system, therefore, one 
might expect the universal scaling function for an isotropic system to be that observed in 
the asymptotic critical region. We assume that the magnetic coupling is weak enough for 
fiuctuations in the vector potential A[r) to be neglected. In fact, it is only in this approxi- 
mation that critical behaviour is to be expected |3^. One barrier to this calculation is that. 



in the low-field regime, all the the Landau levels must be included, and the eigenfunctions 
are extremely inconvenient to deal with. Here, we exploit an integral representation of the 



order-parameter propagator [^ in which the sum over Landau levels is carried out once 
and for all. 

Our initial attempts to estimate the scaling function C{x) made use of perturbation 
theory, which yields accurate results for the critical exponents. Unfortunately, perturbation 
theory does not yield a well-controlled approximation to the function C{x), because it gives 
a spurious singularity in the neighbourhood of Tc2{B) and we have been unable to cure 
this problem satisfactorily by ad hoc methods. The alternative we adopt here is to study 
a generalized Ginzburg-Landau- Wilson superconductor having an order parameter with N 
complex components and to make use of the expansion in powers of 1/A^. The scaling 
function must be extracted by means of a renormalization-group analysis. Curiously, we 
have not found in the literature a formulation of the renormalization group that is well 
adapted to the use of the 1/N expansion as a means of approximation. We have therefore 



developed a suitable formulation, which is presented in detail for the zero-field case in p3 
The extension of the 1/A^ formalism and the renormalization-group analysis to the case 
of a nonzero magnetic field are summarized in sections II - IV below. The calculation of 
the scaling function requires a numerical means of estimating several cumbersome integrals, 
and the techniques we have devised for doing this are described in section V. As has long 
been known in connection with the estimation of, the convergence of the 1/N expansion 
is very poor. The next-to-leading order calculations reported here do not yield meaningful 
results for small values of A^ (in particular for the physically relevant value A^ = 1), but 
for larger values we obtain scaling functions which are qualitatively similar to that observed 
experimentally. These results are presented and discussed in section VI. 



II. THE 1/A^ EXPANSION 



We consider the Ginzburg-Landau- Wilson theory for an isotropic superconductor with 
A^ complex order-parameter components (piir) in a fixed, uniform magnetic field of strength 
Bq. It is defined by the effective reduced Hamiltonian density 



N \ / ^ 

^ = E {l(V - iA{r)) 0,(r)|^ + to |0.(r)|^} + j^ E \Ur 



2=1 



(3) 



\i=l 



where to is taken to be linear in temperature (to oc T — Tq), and the coupling strength 
Aq to be temperature-independent. A convenient choice of gauge for the magnetic vector 
potential is A{r) = Bq{0,x, 0), corresponding to a uniform field in the z direction, and we 
have absorbed the charge of a Cooper pair into the magnitudes of A and Bq. As explained 
in [Q, a standard integral transformation of the Hubbard-Stratonovich type allows us to 
express the partition function as a functional integral over an auxiliary field \E', 



N 



(4) 
(5) 



Z = ^/ f llV(j)iV(l)* exp - fd^rH 

i=l 

■l?^exp[-iVif,ff(^)] , 

where the effective Hamiltonian is 

i/eff(^) = /rfVA^^^2(r)-Tr,,,,lnA(r,r';^) (6) 

and M is an irrelevant normalization constant. The propagator A(r, r', \1') is the solution of 
-V^ + 2iBoxdy + Blx"^ + to + i^>{r)\ A(r, r'; ^) = 5{r - r') . (7) 



As in 1^, it is helpful to formulate the 1/A^ expansion in terms of a self-energy to, which 
can be defined as follows. The full 2-point function 



can be expressed as 

Xkz,(T,n\ jXkz,cr,n\ J 



NH,ff{^) 



G''\ry) = j-^B,Y. 



dkzda 



(8) 



(9) 



where the Xk^ o-n(^) ^^^ eigenf unctions of the differential operator in (0), whose eigenvalues 
are the Landau levels E{n, k^) = k1 + (2n + 1)Bq + to, and we define 

to = r(2)(o,o). (10) 

In the limit Bq -^ 0, this agrees with the definition adopted in |]53|. Because the Lan- 
dau eigenfunctions are extremely inconvenient to deal with, we shall exploit the fact that 
^-t{x+x ){y-y )Bo/2q{2)^^^ ^/^ jg g^ translatioually invariant function to write 



d^k 



Qm(^j., r') = e*(^+^')fo-S/')Bo/2 f " ^ g^fc■(r-r')(^(2)(^y^^) ^ 



Using the eigenf unctions given in |EB[, we find 



[r(^)(0,A:. 



(nBo)-^ j dKdkye-^^'+^l^l^''G^^\k) . 



(11) 



(12) 



Owing to the factors of N muhiplying ifcff(^) in (|§) and (|^), the 1/A^ expansion is gen- 
erated by the method of steepest descent. We expand \l/ about the the position-independent 
saddle point by writing 



^(r) = i{to + Bo- to + N~'5) + {2Ny'/^ij{r) 



(13) 



where S is defined by the condition (ipir)) = 0. The propagator A(r, r'; \E') can be expanded 
as A(r, r'; \1/) = A(r, r') + 0(A^~^/^), where the leading term is the solution of 



-V' + 2iBoxdy + B^ + to - Bo A(r, r') = (5(r - r') . 



(14) 



In real space, the diagrammatic expansion is identical to that explained in I^S], to which we 
refer the reader for details, except that the propagators are modified by the presence of the 
magnetic field. The (p propagator A(r, r') is given by 



A(r.r') 



J{x+x'){y-y')Bo/2 



d^k 

(2^ 



^ik-{r—r') 



A(fc) 



(15) 



where, as obtained in [^, A(A;) has the integral representation 



A(fc) 



du (coshi^ow) ^exp —{kl + to — Bo)u — {k^ + ky)r{u) 



(16) 



with t{u) = Bq tanhi?oM. The ijj propagator D{r — r') is translationally invariant. Its 
inverse is 



D-\r - r') = XQ^5{r - r') + -A(r, r')A(r', r) 



(17) 



and its Fourier transform is 



D{k)= fA^i + n(fc) 



where 



1 r d^k' 



An J 2 Jo smh.Bo{u + u' 



X exp 



uu 



T[U)T(U 



U + u' 



T[u) + T[U 



(18) 



^{ki + k'y)-{u + u'){to-Bo) 



(19) 



To make use of this expansion, we need to determine the counterterm 5 introduced in 
(|13|) and the relation between the self-energy tg and the variables to, Aq and Bq with which 
we started. Consider first the expansion for the 2-point function G^'^\r,r') shown in figure 
1(a). The first term is just A{r,r'), which contains the exact self-energy to, and satisfies 
(p!0|) by itself. Thus, the counterterm 6 is required to cancel the one-loop contribution at 
n = kz = and we find 



6 



2nBo 
1 f d^k 
2 



(27r) 



dkkdkye-^^'+^y^^'''' 
A{k)D{k), 



d^k' 
(2^ 



A(k + k')D(k'] 



fcz=0 



where 



A(fc) 



du exp —{kl + to)u — {2Bo 



-i fl _ g-2Bo« 



\kx + ky) 



(20) 



(21) 



The requirement that {ip{r)) = yields a constraint equation, which implicitly determines 
to. Figure 1(b) shows the expansion of (ipif)) to order 1/A^; the function / is the coefficient 



of ^l>{r) in Hes, as given in |]33|. From this we obtain 



to = $o(to, Ao, Bo) = to - 5o - y A + N'^ 



Ar 



A - 5[1 + XoUiO)] 



(22) 



where 



^(^0,-^0, -Bo) 
^(^0,-^0, -Bo) 



^'^ Ar^^ /"M'^'r /•- expHto-BoH 



d^k 

^2^ 



As{k)D{k). 



(23) 
(24) 



The function A3(fc) corresponds to the loop of three propagators in figure 1(b), and is 
defined by 



Asik) = /'rfVdV'e''=("'-"")A(r,rOA(r',r")A(r",r). 



(25) 



Straightforward but tedious algebra suffices to show that it is independent of r and is given 
by 



As{k) 



dU{k) 
dio 



(26) 



Our aim is to investigate the scaling properties of the specific heat. Within the Ginzburg- 
Landau- Wilson approximation, the specific heat per unit volume per order-parameter com- 
ponent is given by 



C 



1 dHnZ 

2NV dtl 



{2N)-'Jd'rj:{mO)\'\Ur)W, 



(27) 



where V = J d^r is the volume and (. . .)c denotes the connected correlation function. This 
correlation function can be obtained directly as 



but it is not hard to obtain the convenient expression 



NH^ 



c = \,' 



1 - Ao-^I?(0) 



(28) 



(29) 



where 'D(k) = D(k) + 0(A^~^) is the Fourier transform of the 2-point function {ip{r)il}{r')). 
To order 1/A^, this 2-point function is given by the sum of diagrams shown in figure 2, and 
is conveniently expressed in terms of a self-energy 11^ (fc) as 



V{k) = D{k) + N-^D{k)U^{k)D{k) . 
As explained in appendix A, the self-energy at fc = is given by 



IdA dU{0) d 



-A-6U{0) 



dt< 



-um 



(30) 



(31) 



the second expression being convenient for the purpose of renormalization. 



III. RENORMALIZATION 



The scaling behaviour of thermodynamic functions emerges in the usual way from a 
renormalizat ion-group analysis, but in the context of the 1/N expansion this requires a 



nonstandard renormalization scheme, which is developed in detail in |Q for the theory with 
Bq = 0. According to this scheme, renormalized variables t, t, z and B are defined by 



tn 



5n 



fih 



66 






z 



l_iV-i^lnz + 0(N-2) 



yU ^Z 

fx^B. 



AC 



(32) 

(33) 

(34) 
(35) 

In these expressions, toc is the value of to at the zero-field critical point, S^ = (27r)^^ 
is the usual factor arising from angular integrations, and the constants a = l/lGvr and 
b = 1/16 arise from the large- and small-momentum limits of n(fe) when B = 0. As 
usual, n is an arbitrary renormalization scale, with the dimensions of inverse length. The 
magnetic field requires no renormalization; the definition (^51) simply serves to make B 
dimensionless, as are t, t and z. In this scheme, critical behaviour is governed by an infrared- 
stable renormalization-group fixed point at z = 0. The criterion for renormalization is that 
renormalized thermodynamic functions should have finite, non-zero limits as z ^ 0, and we 
have implemented this requirement by a 'minimal subtraction' of the leading singularities 
proportional to In 2;. It is crucial to our analysis that, as in the perturbative renormalization 



of [^ , the presence of a magnetic field introduces no additional divergences beyond those 
encountered at B = and we shall return to this point shortly. 

For our immediate purposes, we need renormalized versions of the constraint equation 
(^) and the specific heat (^91). The various integrals and subintegrals from which these are 
constructed must be re-expressed in terms of the renormalized variables. To this end, it is 
convenient to introduce the dimensionless quantities 



Dnip; a, z, B) = [z + B-^/^Unip; a)]"' 

An{a) = ifi'B)-'/^ [A{f,H,fi'B) - A(0,0) 

An{p;a) = i2^BA{k;fiH,i2^B) 



AR{a,z,B 



fi~^ [A{fiH,^z~\fi'^B) - 46{^H, ^iz~\ n'^B)U{0; fiH, fi'^B) 
6R{a,z,B) = ifi^B)-^/'^ [6{ijH, fiz-\ fi^B) - 6{0, fiz-\0) 
A3r(p; a) = {fi'Bf/'A,{k; fiH, fi'B) , 

with rescaled variables p and a defined by 

p = {^^BY^'^k and a = i/B . 

Subsequently, it will also be helpful to write 

2 2 2 a 2,2 2-2/) 

Pz = P COS U , p^ + Py = p Sin & . 

With this notation, the constraint equation becomes 



t = {z + 2a)-^<^{i,z,B) 



mz,B) 



z{i-B)-^-B'/'A^ 



+N-^ 






-zB^lH^\[^-t-^-B\^z\nz 



2^ 



2 
3' 



+ 0(N- 



The dimensionless, renormalized specific heat Cii{t, z, B) is defined by 



C(to, Ao, Bo) = C(0, Ao, 0) + CiAo-'(to - toc) + Z^Cnit, z, B) , 



where 



Zt 



[z + 2a) 



l_Ar-i^ln2 + 0(N-2) 



(36) 
(37) 
(38) 
(39) 
(40) 
(41) 
(42) 

(43) 

(44) 

(45) 



(46) 



(47) 



(48) 



is the renormalization factor appearing in (33). The dimensionless constant Ci multiplies a 
non-singular term, whose presence in three dimensions was first noted by Abe and Hikamai 



and whose role in our renormalization scheme is discussed in [^. Writing CR(t, z, B) 
z + 2a)^CB.{i,z,B), we find 



8 



Cn{i, z, B) = -D^iO- «, z, B) - -Ci$(t, z, B) 



+N-^ Ei{i,z,B) + E2{i,z,B)+E3{i,z,B) + 0(N 



where 



Ei = -D^B — DK—-\nz-D^B 'a— — — In 

5<5r 5: 

r" 



E2 = DlB-'Un 



E-i 



da 
D^ — —zliaz 



^z-^ - Dn^ \nz-Ci{ z{i - B) - ^B^'^^k \^\nz 



■26 



53 
6 



'R 



36 



(49) 

(50) 
(51) 
(52) 



and Dr and IIr stand for -Dr(0; a, z, B) and nR(0; a) respectively. The integrals Ar(q;, z, B) 
and (5R(a, z, B) defined in (^) and (^) are both singular when 2; — > 0, but each of the 
quantities $(t, z, -B) and Ei{t, z, B) has a finite limit, provided that we choose Ci = —2/6^. 
To verify this assertion is not an entirely trivial matter. In particular, to verify that the 
expression (H6|) for $(t, 2;, B) has a finite limit, it is necessary to show that 



A^{a,z,B) = -B^/^ U^^{a) 



a- 



9ArM\ & 



36 



da 



ln2 + O(zlnz) . 



(53) 



Because the singularities arise from the large-p region of integration, the required cancella- 
tions can be verified by means of large-p expansions of the subintegrals nR(p; a), A3r(p; a) 
and Ar(p; a), which are discussed in appendix B. 

IV. RENORMALIZATION GROUP AND SCALING 

The fact that the unrenormalized theory is independent of the renormalization scale ix 
leads in the standard way to renormalization-group equations for the renormalized quantities 

t(t, z, B) and CR(t, z^ B), which take the form 



■A.)|-(2-,(.-))f|-2BAH--ij 



■'"4-(-"«4--^-^ 



t(t, z,B) = 
CkH z,B) = 



(54) 
(55) 



where a{z) = 2 — 3iy{z) and the remaining functions are those derived in [Q. In contrast to 
perturbative renormalization schemes, the additive renormalizations of the specific heat in 
(^) are independent of fi, so the associated renormalization group equation ( ^5]) is homo- 
geneous. The asymptotic critical behaviour with which we are concerned here is governed 
by the infrared-stable fixed point at z = 0, where /5(0) =0 and the other function reduce to 
the critical exponents 



16 



+ 0(N-2), 



16 



a 



+ 0(N- 



V 



+ 0(N- 



For z = 0, the renormalization-group equations are equivalent to the relations 



(56) 



9 



t{i, 0, B) = i^i^tir^^-^H, 0, r'^B) = £i/'^(2a)-^$(r(2-^)t, o, r^s) (57) 

c^ii, 0, B) = r"/^CR(r (2-^)t, 0, r^B) , (58) 

where i is an arbitrary scaling factor. The functions t{t, 0, B) and CR(t, 0, B) have infrared 
singularities when their first argument t vanishes. On the right-hand sides of (^) and (p8|), 
we exponentiate these singularities into the prefactors by choosing i to satisfy i~^'^~'^H = 1. 
Then, by setting £ = B^^'^L, we find that Cr has the scaling form 

where, with x = tB^^^'^'^, the scaling function is 

C{x) = L-/^(a;)CR(l, 0, L-\x)) , (60) 

the function L{x) being determined by the constraint equation 

2ax = L^/''{x)^{l, 0, L-^{x)) . (61) 

To obtain a numerical estimate of the scaling function C{x), we need an approximate means 
of evaluating the integrals in (^61) and (^) which is discussed in the following section. 

V. NUMERICAL ESTIMATION OF INTEGRALS 

In order to determine the functions Cr(1, 0, L""^) and $(1, 0, L'"^), and hence the scaling 
function C{x), we need to estimate the renormalized counterparts of integrals such as (pOD 
and (P^. This requires analytic approximations to the functions Ilii{p;a), Ar(p; a) and 
^3r{P]<^), which are themselves defined by rather intractable integrals. This section indi- 
cates the methods of approximation we have used, focussing on the example of A3r(p; a), 
which we express in terms of the variables ( ^41) as A3r(p, 6; a). Having used the renormal- 
ization group to replace i with 1 and B with L~^ in (|60|) and (0), we have a = i/B = L^. 

For large values of p, we approximate all of the subintegrals by means of the large- 
momentum expansions developed in appendix B. Numerically, this turns out to be a good 
approximation for p > 6. In particular, this strategy allows us to cancel analytically the 
divergences that arise at the fixed point z = 0. 

For p < 6, an expansion in inverse powers of a is possible when a is large enough. More 
specifically, in the case of A3R, we have an expansion of the form 

A3r(p, e- a) = «-3/2 [fo{q, 9) + a- Vi(g, 0) + «-V2(g, 9) + 0{a-^)\ (62) 

where q = a~^''^p. (This entails, of course, a rescaling of the integration variable in the final 
integral (0).) Using the representation (|B4|) , the change of variables v = w/a, v' = w' /a 
leads to a power series expansion in which each of the remaining integrals can be calculated 
analytically. We find that the functions fi{q,0) are given by 

/o = (87r)-ig, f, = (l6n)-\Q + 8Q') 

/2 = (967r)-i [3Q + lag" + (128 + 4:qh^)Q^ + 96qh^Q^] , (63) 

10 



where Q = {q^ + 4)~^ and s = sinO. In practice, we have used this approximation for 
a > 2.25, where it appears to yield results of satisfactory accuracy. 

For p < 6 and a < 2.25 no systematic expansion in any small parameter will serve our 
purpose. Instead, we have devised an approximation scheme which we again illustrate for 
the example of Asr. The basic strategy is to evaluate the double integral ( |B2[ ) numerically 
for selected values of p, 6 and a and to construct an interpolating function from these nu- 
merical values. To interpolate simultaneously in all three variables is a difficult undertaking, 
however. To simplify it, we introduce a further approximation, which reduces the function of 
three variables to several functions, each depending on only two variables. In the expression 
^, we make the change of integration variables 



u 



pcos 



u 



psm 



(64) 



The integral becomes (again, with the notation s = sin 6 and c = cos ( 

„3/2 



A3r(p, 6; a) 



np, 



1 /•i'/2 j-oo p^ 

^,, ,n/9 / dcj) sin(2(/)) / dp ^— — 
2(47r)3/2 7o 7o smhp 



X exp 



a-l)p- (pV/4)psin^(20) -p^s^T{p 



tanh(p cos^ (p) tanh(p sin^ <; 
tanh(p cos^ (p) + tanh(p sin^ 



(65) 
(66) 



Our approximation scheme is based on the observation that T{p,(f)) ^ (p/4) sin^(20) for 
p — > while T{p, 0) ^ | for p ^ cxd, except at the endpoints = and = 7r/2. We divide 
the region of integration into two parts: region I, where < p < 5(0), and region II, where 
S{4>) < p < oo. The boundary p = S'(0) is determined in a manner to be explained shortly. 
We have A3R = A\^ + Af^, where 



.1 



1 



n3/2 



.7r/2 rSW p- 



X exp 



A" 

^3R 



-p2/2 



rTv/2 roc p 

o/. N^/9 / (i0 sin(20) / dp^- — 
2(47r)3/2 7o Jsw smhp 



a - l)p - (/p/4) sin^(20) + p's'X' 
3/2 



X exp 



2„2x^II 



-{a - l)p - (p^c74)(psin^(20) - 2)+p's'X 



(67) 



(68) 



with 



X\pA) 
X"(p,0) 



(p/4)sin2(20)-r(p, 

h-np,4 



(69) 
(70) 



and we propose to expand the integrands of (|67|) and (^) in powers of X^ and X^^ re- 
spectively. The boundary p = S'(0) is chosen as the locus on which X^ = X^^, namely 
5(0) = 2/sin^(20), so that the two expansions match term by term on the bound- 
ary. With this choice, X^{p, 0) and X^^{p, 0) are always smaller than the boundary value 
X^{(j)) = i — T{S{(f)), 0), which itself has a maximum value of approximately 0.18 at = 
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and = '7r/2. Moreover, the quantities p^s^X"^ are positive, so the expansion of each inte- 
grand converges monotonically. This is not necessarily true of the integrals, but in practice 
we have found that retaining only the first two terms of each expansion yields results that 
are fairly accurate and match smoothly to the large-a and large-p expansions. It will be 
seen that each term in the expansion of /S}^^{p,6]Q.) is of the form {p'^s'^)"'f^{p'^,a), while 
each term in the expansion of Aj^^iiP^ ^5 '^) is of the form e~^ ^'^{p'^s'^)"'fn{p'^c^, «)• Each of 
the functions f^ depends only on two variables. We have obtained interpolations for these 
functions, for tt, = 0, 1, giving final approximations for A^r^ of the form 



A3r(p,^;«) 



A^kfe^;«) 



r Rl,{a 


) 1 


1/6 

1 „2,2 




1 + Rl, 


^ P 





^2,0 («) 



1 + RY^2{a)p^c^ + RtM)P^c'^ 



ELl^2,n(«)F 



^ 2 2 e-^ /'^2!o(«) 

+ p s 



,2n 



1/2 



l + i?^ii(a)pV 



(71) 



(72) 



where the Rf^ are rational approximants obtained from the Thiele interpolation formula. 



I, J 



The form of these interpolating functions (and those for nR(p; a) and Ar(p; a), which we 
do not give explicitly) is chosen so as to give the correct behaviour at large values of p and 
to allow the integrals over 6 to be performed analytically in the final calculations of A-^ and 

VI. RESULTS AND DISCUSSION 



As compared with perturbation theory (the expansion in powers of Aq), the 1/A expan- 
sion has a key advantage when applied to the problem of a superconductor in a magnetic 
field. At the lowest order of perturbation theory, it is not hard to show that the specific 
heat behaves as 



C ~ Bo{to - ioc + Bo 



-3/2 



5o[T-Te2(5o)] 



-3/2 



(73) 



in the neighbourhood of the line T = Tc2{Bq). In terms of the formalism used in this paper, 
this approximation corresponds to taking C ~ 11(0) in ( ^91) and to ~ ^o ~ toc + -Bo in the 
constraint equation (P^). This divergence is entirely spurious. Experimentally, there is no 
sign of it and theoretically it can be removed by means of a self-energy resummation of the 
Hartree variety. The 1/A^ expansion incorporates this resummation in a way which allows a 
renormalization-group analysis of the scaling behaviour to be systematically pursued. 

On the other hand, the 1/A^ expansion has serious drawbacks. Even at the next-to- 
leading order we have used here, calculations are extremely cumbersome, and this is unfor- 
tunate, because the convergence of the expansion is notoriously poor. With the relevant 
value A^ = 1 for the number of complex order-parameter components, the specific heat ex- 
ponent given in (|56|) is a ~ 0.62; compared with the best theoretical and experimental value 
for the XY model axY ~ —0.013, it is about fifty times too large and has the wrong sign! 
For the correlation-length exponent, we have u ~ 0.46 compared with uxy ~ 0.67. 

Using the formalism and numerical methods summarized above, we have obtained esti- 
mates for the specific-heat scaling function C{x) as given in (^). Here too, we find that 
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the convergence is poor; for small values of N, the next-to-leading terms are larger than the 
leading terms. Matters are somewhat improved if the XY exponents z/xy and axv are sub- 
stituted in (^) and ( |6TD for those shown in (^). Here, we present only the best results (as 
judged by their qualitative similarity to experimental data) that we have been able to obtain 
by this strategy. Figure 3 shows the scaling function calculated for values of A^ between 10 
and 20 and, for comparison, figure 4 reproduces the experimental data reported in |^. For 
more negative values of x than those shown in figure 3, the calculated curves diverge rapidly, 
either to large positive values or to large negative values, and our approximations are clearly 
inadequate in this region. The reason for this is not entirely clear to us. One possibility is 
that our neglect of the non-zero order parameter {4>{r)) in the mixed state becomes seriously 
inadequate at temperatures a little below Tc2{Bq). In the lowest-Landau-level approxima- 
tion, it seems to be possible to continue the scaling function of the homogeneous normal 
state to temperatures well below Tc2{Bq), but the same may not be true of the critical-point 
scahng function. 

Since we cannot obtain reliable results for the physically relevant number of order- 
parameter components A^ = 1, a detailed fit of our calculated scaling function to the data 
would have little meaning, and we have not attempted it. For larger values of A^, it is 
clear that the calculated scaling function does reproduce the qualitative features of the ex- 
perimentally determined function in the region where our approximations appear to work. 
Tentatively, at least, it seems reasonable to conclude that the critical-point scaling implied 
by the Ginzburg-Landau- Wilson model is consistent with what is actually observed. A less 
tentative conclusion is that some much better method of approximation than those currently 
available is needed to test this scaling prediction quantitatively. 

Whether the scaling observed in YBCO really corresponds to a regime dominated by 
critical-point fluctuations is another matter. Indeed, the scaling functions exhibited in figure 
3 are qualitatively very similar to the 3-dimensional scaling function of the lowest-Landau- 
level approximation estimated by Tesanovic and Andreev ||5^. It is far from clear, therefore, 



that a quantitatively more reliable estimate of the critical-point scaling function, should it 
be obtainable, would serve to discriminate between the two scaling hypotheses. 

In this paper, we have attempted to estimate the scaling function associated with asymp- 
totic critical behaviour, which is controlled by the renormalization-group fixed point z = 0. 
In principle, the formalism and numerical approximations described here should also facili- 
tate an investigation of the competition between low-field critical-point scaling and high-field 
lowest-Landau-level scaling in the intermediate region where neither type of scaling is exactly 
valid. We plan to address this issue in a future publication. 
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APPENDIX A: CALCULATION OF VERTEX FUNCTIONS 

As explained in detail in [Q, the basic elements of Feynman diagrams in the 1/A^ 



expansion are vertex functions of the form 

A„(ri, . . . , r„) = A(ri, ra) A(r2, rg) . . . A(r„, r^) . (Al) 

These are gauge- invariant functions, and therefore also translationally invariant, but in the 
presence of a magnetic field, the form of the propagator (|I5|) makes them somewhat awkward 
to handle. For the purposes of this paper, the fact that we usually need to integrate over 
one or more of the arguments rj leads to some simplification. Let us write 

A(r, r') = e*(-+-'){j/-j/')Bo/2^(^ _ ^/) ^^2) 

where A(r) is the function whose Fourier transform is given in ([T6|). We find that 

/ £r"A{r, r")A{r", r') = e^(-+-')(s/-s/')Bo/2^^(^ _ ^/^ ^ ^^3) 

with 

A2(r) = / d^r' e'i-'y-y'-)Bo/2^^^ _ r')A{r') . (A4) 

A lengthy, but straightforward calculation shows that the Fourier transform of A.2{r) is 

A.W^-^. (A5) 

dto 

When Bq = 0, this reduces to the familiar fact that d[{k'^ + io)~^]/dio = — (A;^ + to)"^- The 
function A^i^k) defined in (p5|) is equivalent to 



AJk) = f -^ A(k')A2(k + fc') = --i / 7^ A(k')A(k + k') = _S^ . (A6) 

In the same way, we can define A^lr — r') by 

J dh"£r"' A(r, r") A(r", r'") A(r"', r') = e^(-+-'){i/-i/')Bo/2^^^^ _ ^/) ^ (^^^ 



and find that its Fourier transform is 

A3W = -i^ = i^. (A8) 

Consider now the self-energy diagrams shown in figure 2, which are to be evaluated with 
the external wavevector equal to zero. The first one is 
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and the second is 



<^(0) 



(27r)3 (27r)3 
1 r d^k d^k 



D{k)A{k')A3{k + k') 

,.dA2ik + k') 



so we can use f|A6|) to write 



2 7 (27r)3(27r)3 



r(2)/ 



D{k)A{k' 



dto 



(AlO) 



<)(o) + 2n;'=^(o) 



d^k ^^^^OAm 



(27r 



dto 



(All) 



Using the expression ([TsD for -D(fc) and the first expression in (|A6|) for A3(fc), we find that 

dD{k) 



dtn 



-D{kyA,{k) 



d^k 



so the third diagram of figure 2 is 

^f(0) = J-^^A;{k)D{krA,{k) 
f d^k dD{k) ^ ,,, 

and we obtain 



r(i). 



r(2)/ 



r(3)/ 



d r d^k 



U\;>{0) + 2U\;>{0)-U^^>{0) = -— I j,^^D{k)As{k) 



dto J (27r)3 



dAiio,Xo,Bo 
dto 



The final diagram in figure 2 is 



uf{0)=As{0) 



dU{0) 
dto 



(A12) 



(A13) 



(A14) 



(A15) 



APPENDIX B: LARGE-MOMENTUM EXPANSIONS 

To verify that the constraint equation and the specific heat can be correctly renormalized, 
and also to assist the numerical estimation of the renormalized quantities, we require expan- 
sions of the subintegrals nR(p; a), A3R,(p; a) and Ar(p; a). We use the notation indicated 



in (|4^) and the abbreviations s = sin 6* and c = cos 6*. For Ar(p; a), the expansion 



A^^{p,9;a) 



du exp 



-(pV + a)M--(l-e-2>V 



P 



-2 



(a - 2s^)p~-^ + 



a 



(6a + 4)s^ + 12s^ 



P 



0(p 



-8\ 



(Bl) 



follows trivially from the change of variable u = v/p^. For A3r(p; a) we have the expression 
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^3R 



2(47r)3/2 Jo 



{u + u'Y^^exp 



du du ■ 



-{a- l){u + u') 



UU 2 2 '^'^ 2 2 



{u + u') 



7\P ^ 



[T + T 



■kP s 



sinh(-u + u') 



(B2) 



where r = tanh-u and r' = tanh-u'. By virtue of the symmetry of the integrand under 
interchange of u and u', the region of integration < -u' < -u yields exactly half of the 
integral. In this region, we can make the change of variable 



u + u' = V + v' , Auu' /{u + u') = v' 



to obtain 



4(47r)3/2 Jo 



1 roo {v + v')exp —{a — l){v + v') — jv'p^c^ — (j{v,v')p'^s'" 

AsR = ., . . „ ,^ / dv dv' 



t)i/2 sinh(f + v') 



where 



a[v,v 



cosh(f + v') — cosh ( Jv{v + v' 



/2 sinh(f + v') . 



(B3) 



, (B4) 



(B5) 



A further change of variable v' = v" /p facilitates an expansion in powers of j9 , with a 
result of the form 

A3R(p,6';a) =(5o(a)p"^+ Qi{a) + Q2{a)s^ p"^ 



+ 



QM) + QM)s' + Q^{a)s* P"' + 0(p 



The coefficients Qi{a) are 
with 



(B6) 



(B7) 



Qo('i', tt) = 1 5 Qilf, a) = 4(f ^ — cothf + 1 — a) , Q2('i^, a) = 2(cotht; — f ^) 
Q3(t;, a) = 16 [a^ - 2a + 2(cotht; - 1 + a)(cotht; - v'^) 

Q^{v,a) = 4[5 + 3(cotht;-t;~^)(t;"^ - 4cotht; + 2 - 2a) 
Q5{v,a) = 12(cotht;-t;-^)^ 

The function Il^{p;a) satisfies dU^^ip; a) / da = — A3R(p;a), but does not itself have an 
expansion in powers of p~^. At -B = 0, we have the exact result 

U{k; i, 0) = (SttA;)-^ tan-\k/2i^/^) = bk~^ - Aai^^^k-'^ + 0(k-^) , (B8) 

with k = \k\, which implies that IIr{p', a) has the limiting form 

Un{p; a) = bp'^ - Aaa^'^p'^ + 0(p-^) (B9) 
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as a ^ oo. By integrating A3r(p; a) with this boundary condition, we obtain the expansion 

nR,(p; a) = bp~' + AR,(a)p-2 + 0(p-^) , (BIO) 

which is sufficient for our purposes. 

In the constraint equation (|46| ) and the specific heat (|49|) , singularities at z — *> arise 
from the large-p region of integration in integrals of the form 

/ Wf ^^^^' "' ""' ^^^^^' "^ • ^^^^^ 

By restricting the range of integration to |p| > po, where the value of po is immaterial, 
the leading singularities can be extracted by means of the power-series expansions given 
above. Using the expansion D^^{p; a,z,B) = z + B'^^'^bp^^ + B^^/^AR(a)p~^ + 0(p~^), we 
encounter the three divergent integrals 

^'^ ' -5v4ln. + ... (B12) 



(27r)3 {z + B-y^p-^)p^ b 



dp 1 ni/2'S'; 



(27r)3 {z + B-y^bp-^Yp* b 

d^p 1 .,^Ss 



By^^z-' + ... (B13) 



fi=^/2-|ln^ + ... , (B14) 



(27r)3 {z + B-^/%p-yp^ 63 

where the ellipsis represent less singular terms. These results, together with straightforward, 
though tedious, manipulations of the integrals (|B7| ) suffice to verify that the functions $, 
El and E2 have finite limits and, in (HBf), that Z(5r, = — Ar5'3/26 + . . .. 
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FIGURES 

FIG. 1. Diagrammatic representation of (a) the order-parameter 2-point function and (6) the 
expectation value {tp{r)) at next-to leading order. 

FIG. 2. Diagrammatic representation of the 2-point function for the auxiliary field ^ at 
next-to-leading order. 



FIG. 3. Numerical results for the specific-heat scaling function C{x) for several values of A^. 
FIG. 4. Experimental data for the specific-heat scaling function as reported in Ref. [4]. 
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